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We investigated the SU(2) Einstein- Yang- Mills system on a time-dependent non-diagonal cylin- 
drical symmetric spacetime. From the numerical investigation, wave-like solutions are found, con- 
sistent with the familiar string-like features. They possess an angle deficit which depends on the the 
initial form of the magnetic component of the Yang-Mills field, i.e., the number of times it crosses 
the r-axis. The appearance of a curvature singularity depends on the ratio of the Planck scale M p i 
and Yang-Mills couplings constant g, i.e., a = M p i/g = G an d the initial mass per unit length 
, of the cosmic string. The solitonlike behavior of gravitational and Yang-Mills waves show signifi- 

■ cantly differences from the ones found in the Abelian Einstein-Maxwell system. The stability of the 

' resulting non- Abelian spinning string is analyzed using the multiple-scale method. To first order a 

consistent set of equations is obtained. 

However, a throughout analysis is necessary in order to conclude if the solution is topological 
stable and will have magnetic charge. 
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I. INTRODUCTION 



Taking the Yang-Mills model as the matter part of the equations of Einstein, one finds a rich spectrum of stationary 
regular and non-regular solutions. This in contrast with the vacuum solution and Einstein-Maxwell(EM) solutions, 
where the Schwarzschild respectively the Reissner-Nordstrom(RN) solution are the only static black hole solutions. 
In the Einstein- Yang-Mills (EYM) theory one finds not only the embedding of the Abelian RN black hole, but in 
addition there are genuine non-Abelian coloured black holes. Their co-existence gives rise to the violation of the 
no-hair conjecture, since they carry the same (magnetic) charge. For an overview, see Volkov and Gal'tsov [0. 
After the discovery by Bartnik and McKinnon(BMK) H of these non-Abelian particle- like and later the non-RN 



(hairy) black hole solutions of the SU(2) EYM theory jjj, throughout instability investigations where done by several 
authors Q It turns out that under spherically symmetric perturbations the BMK solution is unstable. The nth 
Ch . BMK solution has in fact 2n unstable modes, comparable with the flat space-time 'sphaleron' solution. In the static 
case, it was recently conjectured || that there are hairy black holes in the SU(N) EYM-theories with topological 
instabilities. In order to analyze the stability of the solutions, one usually linearized the field equations (^Q. We 
conjectured that the conventional linear analysis is in fact inadequate to be applied to the situation where a singularity 
is forme d ||10||. One better can apply the so-called multiple-scale or two-timing method, developed by Choquet- 
Bruhat |13|. This method is particularly useful for constructing uniformly valid approximations to solutions of 
perturbation problems [Q. At the threshold of black hole formation it is found numerically Jl5| that there is 
a critical parameter whose value separates solutions containing black holes from those which do not. For a critical 
value one observes self-similarity: it produces itself (echoes) on progressively finer scales (Choptuik-scaling) . It is 
evident that the collapsing ball of field energy will produce gravitational waves, which will be coupled to the YM-field 
perturbations. Due to the accumulation of echoes, the curvature diverges. Close related to this critical behavior is 
the comparable irregularities which can be found in the Einstein-Skyrme(ES) model and the Einstein- Yang-Mills- 
Higgs(EYMH) model. In the latter model, space-like hypersurfaces develop for a critical value of a = M\y /gM p i two 
distinct regions separated by a long throat, with M\y the YM scale and M p i = 1/ V 4arG. 

The solutions mentioned above need not be spherically symmetric. Since the pioneering work of Stachel |J], we 
know that stationary cylindrical symmetric rotating models possess intriguing features. First, the evolution of the 
'classical' cylindrical gravitational wave-solution in the vacuum situation can be described in terms of two effects: 
the usual cylindrical reflections and a rotation of the polarization vector of the waves, an effect comparable with the 
Faraday rotation [ pp[ . Secondly, the time-dependent solutions can be related to the cosmic-string solution interacting 
with gravitational waves H . Usually the cosmic string solution, with its famous conical structure ,is found in the 
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U(l)-gauge Einstein-Higgs system. It is remarkable that the cylindrical rotating vacuum solution possesses the same 
conical structure. Moreover, it is asymptotically flat and free of singularities. The physical interpretation can nicely 
be formulated using the Belinsky-Zakharov method of integrating the Einstein equations with the help of the Ernst 
potentials ^2|. Using the C-energy description, one can conclude that the solution represents solitonic gravitational 
waves interacting with a straight-line cosmic string, where the asymptotic angle deficit depends on the C-energy of the 
gravitational waves. When a Maxwell field is incorporated, the physical interpretation of the solutions is not easily 
formulated. It is hard to believe that in the electrovac solution the electro-magnetic waves do not contribute to the 
angle deficit |2j|. Moreover, the mass per unit length of the cosmic string, which is a constant in the Xanthopoulos 
solution, will be determined by the energy-momentum tensor of the non- vacuum situation. In stead of 'electrifying' or 
'magnetifying' the vacuum solutions via the Ernst formulation, one can better consider the coupled system of Einstein 
equations and matter field equations. 

Here we consider the Einstein- Yang-Mills system on a cylindrical symmetric rotating space-time. This situation 
con formally be obtained from the axially symmetric situation by complex substitution of t — * iz,z — > it ^4| . In 
the non-rotating cylindrically symmetric situation, it was numerically found [^6| that the singular behavior at finite 
distance of the core as found in the supermassive situation, will be pushed to infinity for some critical values of the 
gauge coupling constant and one of the YM components. 

We will use a (2+l)+l reduction scheme in order to obtain the suitable parameterization for the Yang-Mills 
potential. In fact one replaces the YM field by a 2D Abelian one-form, a 2D Kaluza-Klein two-form, a 2D matrix- 
valued scalar and a 3D Higgs field [^?J . 

In the Abelian counterpart model it is known that the space-time around (spinning) gauge strings could have, besides 
the usual angle deficit feature, exotic properties, such as the violation of causality [ p8|j29| , time-delay effects |30|] , helical 
structure of time |3lj-|33| and frame-dragging p4| . A class of approximate solutions of the coupled Einstein-scalar- 
gauge field equations on the (r, z)-plane of an axially symmetric space-time was found fjSfl . So a question addressed in 
this paper will be whether the string features will be maintained in the EYM system. In the general (r, z)-dependent 
situation, a nice overview was given by Islam |36| . One can proof via the Ernst formulation that for any solution of the 
stationary Einstein equations one can generate a corresponding solution of the EM equations. In fact, one introduces 
two complex potentials, for which one finds two second-order elliptic differential equations. Starting with the Kerr 
solution, it is possible to use this correspondence to find a solution of the EM equations leading to the Kerr-Newman 
solution. This 'Ernst-route' fails in the EYM system, simply by the fact that one of the metric components does not 
decouple from the other ones. It is also known that the electrovac staticity theorem |}7| does not generalize to the 
EYM system. In fact, spinning EYM systems must be electrically charged. So another question addressed in this 
paper is whether electrically and magnetically charged solution will exist in the EYM system. Finally, this is not the 
complete story: the investigations can be extended to the model where the interior solution is properly matched onto 
the exterior solution of the string |3q| . This is currently under study by the author. 

The plan of this paper is as follows. In Sec. II we derive the field equation on an cylindrical symmetric space time. 
In Sec. Ill we present the numerical solution. In section IV we present a approximate solution using the multiple-scale 
method and in Sec. V we summarize and analyse our results. 



II. THE FIELD EQUATIONS OF THE EINSTEIN YANG-MILLS SYSTEM 



Consider the Lagrangian of the SU(2) EYM system 

S = J cfix^—g 

with the YM field strength 
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g the gauge coupling constant, G Newton's constant, the gauge potential, and TZ the curvature scalar. The field 
equations then become 



(3) 



V^ va = 0, (4) 

with T the energy-momentum tensor 
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and V the gauge-covariant derivative, 

V a T% = V a J%, + ge abc A b a ^. (6) 
Let us consider the cylindrical space time 

r 2 

ds 2 = -n 2 (dt 2 -dr 2 ) + f(dz + uudifi) 2 + jdip 2 , (7) 

where ft, f and w are functions of t and r. The problem is, which ansatz one needs for the YM field. It is often 
argued that Einsteins gravity on a three dimensional space time shows some similarity with Chern Simons theory 
with the Poincare group being the underlying gauge group. So when trying to find the dimensional reduction for the 
stationary axially symmetric EYM system, one can be guided by the results found by using the (2+1)+ 1 reduction 
scheme suggested by Maeda, Sasaki, Nakamura and Miyama p5fl , and later worked out by Gal'tsov p7| . One writes 
in general 

ds 2 = ~e^{dt + Vl dx 1 ) 2 + e^-^idip + k a dx a f + g ab dx a dx b , (8) 

where Vidx 1 = uj(dip + k a dx a ) + u a dx a , k a an one-form, generating a field strength n a b = d a kb — dbk a . The one-form 
v a give rise to the 2D field strength u> a b = d a Vb — dbV a + wK a {,(a,b=l,2). The YM potential is then written as 

A^dx^ = a. a dx a + $(dt + Vidx 1 ) + &(d<p + k a dx a ) (9) 

Here a is the dynamical part of the YM field, parameterized by a 2D complex Abelian one-form, <I> a 3D Higgs field 
and ^ a matrix-valued scalar. In our simplified case of metric (7), we have v\ = —u>ki,k 2 = v 2 = and the YM 
potential can be as 

A,j, = $(dz + ud(f) + ^{dip - kidt) + a a dx a (10) 
The YM-part of the Lagrangian can then be written as 

2/, 



Cym = t 



{fabf ab + -^(Da* + $d a io - a' a )(D a V + $>d a u - a' a )} 
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j{» a *D"* + 4(*'+9[*,*]) 2 }], (11) 



with 

D a * = da i & + g[a a ,*}-k a y', (12) 

f a b the 2D field strength 

fab = d a a b - d b a a + g[a a , a b ] + a a k b - a b k a + &u ab + ^k ab , (13) 

and a prime denoting the partial derivative with respect to ip. 
Let us consider the following parameterization 

ax = A 2 (t,r)T v , a 2 = Ax{t,r)r v , * = Wx(t,r)r r + (W 2 (t,r) - 1)t z , 

* = $i(t,r)r r + $ 2 (t,r)r z , (14) 

with Tj the usual set of orthonormal vectors. From the condition T u — T rr — 0, we obtain 

d t A 1= d r A 2 , W 2 ='^± + 9 —. (15) 

*i g 

From the Einstein equations and the YM equations we obtain the set of partial differential equations for 
fi, /, cj,$i,$2 an d W\. It follows from a combination of the YM equations, that A% and A 2 can be expressed in 
$1 and <& 2 : 
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For the gauge 9*^2 = 9 r Ai and hence d 2 Ai — d 2 Ai = 0, we obtain from the Einstein equations Eq.(3) 
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(dtf) 2 - {d r ff 



Orf = 0, 



8nG 



d 2 tf-d 2 r f+- f [{d r ff-(d t ff 



~drf = 



r 2 $ 4 



(<i>4 r 2 _ /2$2 w 2 )((5r$i) 2 _ (9t$i)2) + ($ 4 r 2 _ ^^JjJ^j^a _ (^$ 2 ) 2 ) 

+/ 2 *?( ( &? + <f> 2 )((W 2 - (w 2 + $?((a^) 2 - (9^) 2 )) 

+ 2/ 2 $3($2 + $2^^^^ _ 9^9^) + 2f 2 <f>j<I> 1 W(d r Wd r <5> 1 - 5 t W5 t $i) 
+2/ 2 $3$ 2 VF(a t $ 2 5 t w - d r <S> 2 d r u) + 2f 2 <S> 2 1 <S> 2 W(d t Wd t <S> 2 - d r Wd r $ 2 ) 
+2/ 2 $ 1 $ 2 W /2 (<9 r $ 2 <9 r $i - dt^ 2 dt^i) + 2f 2 <5> 2 2 <$> 2 W{d r <$> 1 d r u - d t $i<9 t w) 

+ 3 $ 2 ($ 2 r 2 - / 2 vK 2 )(2A 1 ($ 2 a r $ 1 - $ia r $ 2 ) + 2A 2 ($xa t $ 2 - $ 2 a 4 $!) + , 9 ($ 2 + $ 2 )(a 2 - a 2 )^ 



d 2 CJ - d 2 CJ + -Ad t fd t OJ - d r fd r u) + -d r CJ = [^$1 ((3r$ 2 ) 2 - (ft$ 2 ) 2 ) 

J f J ^ i 

+<f> 2 W(d t <S> 2 d t <f>i - a r $ 2 9 r $i) + <5> 2 <f>l(d r <I> 2 d r Lu - d t <S> 2 d t oj) + $ 2 $i(<9 r $ 2 ,W - d t <f> 2 d t W) 

+<S>\{d r <5> 1 d r W - dt^xdtW) + $?(<9 r $i<9 r w - d t ^ x d t J) 

+g 2 <S> 1 W($ 2 1 + $1){A\ - A 2 ) + 2g$ 2 W(A 2 d t <i> 2 - ^i9 r $ 2 ) + 2g$ 1 <S> 2 W(A 1 d r <S> 1 - A 2 d t $ 1 ) 
and from the YM equations Eq. (4) 

d 2 ^ - o 2 ^ - + ^i^d r Lu) 2 - (9 t c) 2 )+i(a r /a r $ 1 - d t fd t ^) 

f 2 qf 2 $ 2 W 
+ i-(d r ud r W - d t ujd t W) + y ' \ (A 1 8 r L0 - A 2 d t u) 

+ i^(A 1 d r f - A 2 dJ) + 2g(A 2 d t <P 2 - A 1 d r * 2 ) - 9 -^± + g 2 ^ 2 - A%) = 



d 2 ^2 - d 2 $ 2 - + ^((drL0) 2 - (Otcu) 2 ) + j(d r fd r $ 2 - dtfdt*2) + ^(drUjdrW - d t Cjd t W) 

f 2 $ 2 W f 2 W 
+ - r2 ^ 2 (d t wdt$i - d r ujd r $i) + -^-{dr^dr^2 - d t ujdt^ 2 ) + 2g(A 1 d r $i - A 2 d t $i) 

+ ^(A 2 d t f - A x d r f) + 9 -^-{A 2 d t u J - A.drco) + g 2 $2(A 2 - A 2 ) + 9 -^- = 
/ r A r 



d 2 t W X - d 2 r W Y + -0 r W + -Adtfd t W - OrfdrW) + ^-(drfdrUJ - d t fd t uo) 
r f f 

+d t ud t § 1 - d r uud r <S> l + ^p-(A 2 d t W - A x d r W) + g<P 2 (A 2 d t uj - A^cu) 



, n , 2 | gA^ 2 W 16ttG 



+g 2 W(A 2 - A 2 ) + 



w((dt<s>2) 2 - (a r $ 2 



+<S> 1 (d t <I> 1 d t W - d^xdrW) + $1(0^8^1 - drud^x) 
+^ 2 (dt^2d t W - d r ^ 2 drW) + -^(9 r $ 2 9 r $i - 9 t $ 2 9 t $i) + $i$ 2 (a t $ 2 a t w - d r ^ 2 d r uj) 



+2gW<S> 1 (A 1 d r <P 2 - A 2 d t <5>2) + 2gW<$>2(A 2 d t <S> 1 - AtfrQi) + <? 2 $?W($? + $ 2 )(^ 2 - A 2 ) 
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III. NUMERICAL SOLUTION 



In the vacuum situation, the behavior of the cylindrical wave solution is well established [ p9|j4C[| . One can have 
reflection of incoming into outgoing waves with the same polarization. Further, there is a non-linear effect describing 
the interaction between the two different polarization + and x modes. If an outgoing cylindrical wave is linear 
polarized, its polarization vector rotates as it propagates. This effect is often called 'Faraday' rotation and is given 
by 

A B 

tan. 6 a = tan#s = — — (23) 

A + B + 



with 



A+, x = -Adtf ± d r f), B+. x = f -{d t UJ ± drtu) (24) 

r 



/ 

A wave emitted from the axis of symmetry behaves like a shock front and is given by 



/ = w(a 2 + l)/(a 2 + w 2 ), uj = ra(l-w 2 )/w(a 2 + l), Q 2 f = b^/y/P-r 2 , (25) 

with w = {t—\Jt 2 — r 2 ) jr and a and b arbitrary constants. One can generate a two-soliton solution and it is found that 
near the axis log(f2 2 /), which also describes the 'C-energy', does not approach zero. So some matter field must act as 
a possible source of the two-soliton field. This curious conclusion can also be formulated in context of the cosmic string 
with mass unit length m = p/(2p— 2), with p a constant originating from the static solution / = r 2p ,n = ar p ~~ p , 
which on his turn is related to the angle deficit of the string, given by 

5<j) = 2tt(1 (26) 



When a 'Rosen'-pulse In / = J Q r f(f3)dp/^(t - (3) 2 — r 2 is emitted, log(r2 2 /) changes by an amount — k 2 , where k 
depends solely on /(/?). So after the pulse, one has a different string, which let Marder Q to conclude that the 
mass per unit length has decreased by an amount ~ k 2 . Because log(f2 2 /) is related to the 'C'-energy, this confirms 
that the waves carry away energy and are evidently generated by a source on the axis. So a continuous change in the 
parameter p may produce suddenly and inexplicable change in physical meaning. 

Further, it was emphasized in the vacuum case that there is a reflection of the outgoing wave into an ingoing wave 
and visa versa. 

From the Einstein-Maxwell case |2l],|2j| it is known that there exist stable open cosmic strings coupled with gravity 
and electromagnetism. When gravitational and electromagnetic waves are generated, the angle deficit is increased 
asymptotically and near the axis independent of the coupling of the Maxwell field to gravity. Further, it was empha- 
sized that the electromagnetic waves do contribute to the rotation of the spacetime by effecting the 'dragging' of the 
Killing field J^, but not to the angle deficit. 

In our EYM case we solved Eq.(17)-(22) for the initial values 

n(0,r) = C ir p2 ~ p , /(0,r) = r 2p , $ x (0, r) = 1 - , $ 2 (0, r) = e~ C3r - 1 (27) 

For Wi we will take as initial value a node number 1 or node number 2 function, i.e., Wi(0,r) = — tanh(r — 1), or 
Wi(0, r) = 1 — er CiT . Here are C, some constants. 

We used a 100x100 grid and applied the method of lines with cubic Hermite polynomials. The roundoff error 
remained below 0.001. In figure 1 we took w(0,r) = and for W% the node number 2 function. We plotted 

respectively: ft, f, u>, $i, W\, $2, the angle deficit, the gravitational puis wave fl 2 f, g vv , and Wi- In figure 2 we took 
a different initial f2 (p=0.1 in stead of 0.25). We obtain a solitonlike solution. The behavior depends on the initial 
mass per unit length of the string, just as in the vacuum case. In figure 3 we took the node number 1 function for W\. 
We observe that the angle deficit changes significantly as well as the gravitational outgoing pulse. So the Yang-Mills 
waves has an impact on the gravitational waves. In figure 4 we took p=0.1 in stead of p=0.25. We see that the 
behavior of the solution changes significantly: the angle deficit close to the r-axis increases, while the gravitational 
puis 'hangs' close to the r-axes. 

In figure 5 we plotted a long-time run for a different initial oj: an outgoing wave is reflected into an ingoing wave. 

An important question is if the solution admits gauge-charges. The gauge-charges Qe, Qm (electric and magnetic) 
are given by 
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(28) 



with F the YM field strength, *F the dual tensor and the integration over a two-sphere at spatial infinity. In our 
situation, we obtain for the charges 

Qe = J r{ d t {W 1 +u)$> 1 )+gA 2 (W 2 -l + u)$ 2 )+A 2 r r + d t (W 2 + u<$> 2 ) - gA 2 {W x + w*i)] r z }dr (29) 



Qm = J ^{\^{d r ^ 1 +gA^ 2 - f -^{^ l d r ^ + d r W 1 +A 1 +gA 1 {W 2 -l))T r ] 
+ [j(d r ^>2 - gA^) - f -^{d r W 2 + $ 2 d r w - gAiWj] T z )}dr, 
where we used Eq. (15). So in general the solution will possess gauge charges. 



(30) 



IV. AN APPROXIMATE WAVE SOLUTION TO FIRST ORDER 



In the appendix we presented the multiple-scale method for the EYM system. On the metric Eq. (7) we have for the 
tetrad component: = (—1, 1,0,0) and = qj-(1, 1,0,0), which implies that the divergence of the null congruence 
becomes 



1 



(31) 



as it should be. From Eq. (A10) we obtain = -Bf, or, in the notation of Eq. (14), A^ ] = -A { 2 \ From Eq. (A8) 
we obtain 



U l (h . i, \ h h h h U g 22 ^22 + g 33 h 33 

HOI — — "-00 + 111 J, "42 — —llQ 2 , ft-03 — — "13, "23 — 2g^ ' 



First of all we have just as Eq. (15), 



$ 2 VFi 1 
d t A x = d r A 2l W 2 = -i-^ + 1 - - 
$i 9 



(32) 



(33) 



For simplicity we will take $2 = and try to find an approximate wave solution to second order. From the propagation 
equations Eq. (A14) and (A18) we obtain 



+ d v w[ 1] =o, 



(34) 



d v W 2 {1) = 0. 



So we will take ci/ 1 ) and W 2 ^ independent of ip. Further, 



d t A^+d r A\ L) =AY' +a\- 



(i) _ v(2) , v( 2 ) gfti 1 



WW 2 (1) , 



(35) 



(36) 



<M 1} f (1) - - - 

+r + ^(dt*l + W + -±j-{dtf + d r f) 

2 i ^ p<jW 



2r 2 



(37) 



(38) 
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■-{dtu, + d r w)(^ + - — (a t $i + + 

l j 2 Ir 



(d t oj + d r tu)($[ 1] + - + + - ffWf + A 2 ) (39) 



o/(i) _ _ 6(i) _ _ o/(i) 1 

■^jridtf + d r f) + —(d t n + d r ti) + -L-- + -^(k 00 + fcn + 2* l) 



8irGrtf 



(i) 



4/2 



Arf AVL 



(5 t $i + + + wfV&Wx + d r Wi + $i(5 t tD + <9 r u>)) - ffWfViOAi + A 2 ) (40) 



2r 



3r 2 $ 



(41) 



and 



+ wf 3 )((9 t $i + ft$i) + <S>[ 1 \d t W 1 +d r Wx) -g$iW^ \A 1 + A 2 



(42) 



These perturbation equations are linear differential equations of first order in h^u and B^ and contain in general 
the term k^ u and C° In our special case, we see appear in Eq. (38) a second order term k^ v and in Eq. (34) a term 

"(2) . " (2) 

A\ . The remaining perturbation equations will also contain second order term like W-y . 

For the background variables $i, W\ and W 2 we use Eq.(Al3). The resulting equations are similar to the equations 
found in section 2 for $2 = 0. From Eq. (A17) we obtain equations for the background variables f2, /, and ui with 
back-reaction terms, for example, 



o?n - d*n = --d r n 1 

r 



n(d r f) 2 nd r f f 2 Vt{d r Cjf 



n 



2/2 



fr 



2r 2 



2r/ 2 r 2 



[/ 4 (u;W) 2 +r 2 (/W) a ]de + 



+ 



4nG£lf 



f r 2 T 



J 2 {w[ 1] + d^W) 2 + r 2 ($i x) ) 2 + f 2 (W^) 2 



J 



(43) 



So we see that a first order HF gravitational and Yang-Mills wave emitted from the string, will change the metric 
background component Q, just as in the vacuum situation for pure gravitational waves. Moreover, there will be a 
distortion of the shape of the wave during its propagation, which does not occur in the pure Einstein case. In general, 
this is due to the fact that second order terms appear in the perturbation equations of and B®. Now let us 

consider, for example, the equation for W$ i , Eq. (36). For A\ = —A 2 , we can solve W% f° r / = e 2p - F° r example 
for p = 1/4 we have 



W 2 {1) = VrB(t - r)e 



-2/3T 1 



(44) 



which approaches to zero for large r. Here is B(t — r) an arbitrary function of (t — r). Together with Eq.(31) we then 
have to first order 



W 2 {1) =l-l/g+ V^B(t - r)t 



-2/3T 1 



(45) 



which is bounded for large r. One can proceed in the same way for the other first order perturbations by imposing 
some trial solutions for the background variables (!(]] 
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V. CONCLUSIONS 



Non-Abelian gravity-coupled solitons and black holes could have played an important role in gauge theories of 
elementary particles in the early stage of the universe. Besides the well-known spherically symmetric solutions, 
axially symmetric solutions are of interest due to the fact that they admit cosmic strings, which could have played 
the role of 'seeds' for the large-scale structure of clusters in the universe. It is quite evident that the solution found 
here is a radiating one and there is an interaction of the solitonlike gravitational waves and the Yang-Mills waves 
with the rotating string. This can be seen from figure 5 where the angle deficit and Q is significantly changed by 
the pulses. This can also be seen from Eq. (41), where the term on the right-hand side represents the back- reaction 
of the waves on 12. The wave-like solutions found in this model is significantly different from the electrovac solution 
found in the Einstein-Maxwell model. The intervening electro-magnetic waves did not contribute to the angle deficit 
asymptotically, while in our model the magnetic components W% has an impact on the angle deficit for large r (compare 
figures 1 and 3). In particularly the number of times the Wi crosses the r-axis influences the gravitational wave puis. 
Further, for a different mass per unit length of the string, a different behavior of the solution is obtained (compare 
figures 3 and 4). 

In order to conclude if the solution is topological stable, we must solve the perturbation equations to higher order. 
Moreover, a combination of the results of a throughout numerical investigation with the multiple-scale method pushed 
to higher orders, will lead to further understanding of the rich structure of the EYM system on a rotating space time. 
Specially the question if a electrically and magnetically charged solution ('dyon') of the rotating cosmic string exists. 
Further, bursts of high-frequency gravitational waves might be detectable by the planned gravitational wave detectors 
LIGO/VIRGO and LISA, even for strings at GUT scale |42). 

The issues mentioned above are currently under study by the author. 

APPENDIX A: THE EYM EQUATIONS IN THE MULTIPLE SCALE FORMULATION 

Let us consider a manifold M. with two different scales, i.e., a mapping from Ai x 1Z into the space of metrics on 
M: 

(x,O=*0(x,0, xeAUeft. (Al) 

We set £ = wll(x), with II a phase function on M. of dimension of length and lo a large parameter of dimension 
(length)' 1 . 

The multiple-scale method assumes that the metric g^ v and the YM-potentials A a can be written as 

9v» = g»u + - VO^O + -ofc/wO^jO + •••> ( A2 ) 



A« =A% + -m(x'U) + ^C;(x°;0 + (A3) 

where £ = uj\l(x a ) and II a phase function. The parameter u> measures the ratio of the fast scale to the slow one. 
The rapid variation only occur in the direction of the vector l a = Jp.. For a function ^(x" 7 ; £) one has 

where = ^§s-\ifixed and i> = jfjr\ x ° fixed- We consider here the case where the magnitude of the perturbation 
hfj, v of the metric with respect to the background g^ is of order w _1 (for different possibilities, such as for example 
with the leading term in the metric expansion of order lo~ 2 , see p"3[|). Further, we will consider here the hypersurfaces 
II(a; cr ) = est as wave fronts for the background metric g, so (Eikonal equation) 

l a h9 af} = 0. (A5) 

Substituting the expansions of the field variables into the equations and collecting terms of equal orders of lo, one 
obtains propagation equations for B®, C™, h^ and and 'back-reaction' equations for g^ v and A*. It will be clear 
from the propagation equation that there will be a coupling between the high-frequency gravitational field and the 
high-frequency behavior of A a ^ when the singularity will be approached. 



8 



First we substitute the expansions Eq. (A2) and (A3) into the YM-equation Eq. (4). We obtain for the order uj 
equation 



/ / pa 7 7 pti 



= 0. 



(A6) 



For the order adequation we obtain 

->A T-a ^ A 
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^a.J~ fiV + ^cefj^uX ^ai/^A + (jn,a (Fail, "I" ^a/*)^Ay ^^{^v,ot av '""m/)'* 

-h" a l a {l ll B*-l v Bl) = 0, 



witn T ^ = ^g^QuKa - hh^u - lah^u). 



Substituting the expansions into the Einstein equations Eq. (3), we obtain for the order ui equation 

— 1 T a — 1 T a — n 

For the u°-equation we obtain 

R»v + R$ = -8^G{.9 A ^ A + l^Bl - hBDiKp + l v B% - IpBi) 



(A7) 



(A8) 



g^g^iKp + l a B$ - hK){T a aT + i a B a T - i t b%)}, 



(A9) 



with i? M „ the background Ricci tensor and R-ffj an expression in k,hh,ti 2 and h (sec We can simplify the 

equations considerable. If we use Eq. (A5), then Eq. (A6) becomes 



so for periodic we have 



PB% = 0, 



FBI = 0. 



(A10) 
(All) 



Using Eq. (A8) and Eq. (A10) we obtain from Eq. (A7) 



^ 1 - ^au^X + B^Valf! - B^V a l u + Z^V a B" - l v V a B^ + l a (B"^ - Bp >v ) 



+U a h» a B°;-U»C°;+ge abc g 



abc-^ia 



0, 



(A12) 



with V the covariant derivative with respect to the background metric g^ v . Integrating this equation with respect to 
£ yields 



V»J*, = 0. 



(A13) 



(A14) 



Substituting back this equation into Eq. (A12) we obtain the propagation equation for the YM-field 

W(l^ - l v BD + nV^Bl - V V B$ - g^L^x + l a l»h x a B a x 
-U^C; + ge abc r a (2l a AlB c u - l v A h a B^) = 0. 

Multiplying this propagation equation with B va , we obtain the 'conservation'-equation 

\7 a {l a B a v B» a ) = l^h x B va T* x - 2ge abc l^B c v B lJa A b ^ (A15) 

so B®B va is not conserved, unless the right-hand side is zero. Multiplying Eq. (A14) with V , we obtain 

r(l x K-l»h x u )T« x = Q. (A16) 







For the choice l"h ViL = and hence from Eq. (A8), h = 0, the equations simplify again. 
Integrating Eq. (A9) with respect to £, we obtain 

1 



1 ;2^ 8vrG, , 



B aa B a a dt-%TtG% v , 



(A17) 



with r the period of high-frequency components. In the right-hand side we have the back-reaction term due to the high- 
frequency gravitational perturbation, the high-frequency YM-perturbations and the background energy- momentum 
term T. Substitution back Eq. (A17) into Eq. (A9), we obtain the perturbation equations for 



1 



4r 



1 ,■ . 



8ttG 



B aa B^d£ 
1 



BttG 



{l^l u B Xa Bl 



+F\{l v B Xa - l x B a v ) + T a vX (l^ Xa - l X B") + -g^F a ^l a B^ - l?B aa ) 



(A18) 
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Figure 1. Plot of the metric- and YM components for the initial values: £1 = Ar p P (A = 3,p = 0.25), / 
= ^TT $1 = 1 _ e_ °' 8r ' w = 0.8 - e-( r - 2 ) 2 ,$ 2 = e- 32r - l.G=4.5 g=1.25. 

Figure 2. As figure 1, with p=0.1, G=0.6 

Figure 3. As figure 2, for the initial W\: W\ = — tanh(r — 1) 

Figure 4. As figure 3, for a smaller value of the mass per unit length of the string: p=0.09 in stead of 0.25. 

Figure 5. A long time run for $! = 1 — e~°- 5r 
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